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Abelian Chern-Simons-Maxwell theory can emerge from the bosonisation of the 2 + 1-dimensional 
Thirring model that describes interacting Dirac fermions. Here we show how the Thirring model 
manifests itself in the low energy limit of a two-dimensional tight binding model of spinless fermions. 
To establish that we employ a modification of Haldane's model, where the "doubling" of fermions 
is rectified by adiabatic elimination. Subsequently, fermionic interactions are introduced that lead 
to the analytically tractable Thirring model. By local density measurements of the lattice fermions 
we can establish that for specific values of the couplings the model exhibits the confining 2 + 1- 
dimensional QED phase or a topological ordered phased that corresponds to the Chern-Simons 
theory. The implementation of the model as well as the measurement protocol are accessible with 
current technology of cold atoms in optical lattices. 

PACS numbers: 11.15.Yc, 71.10.Fd 



Introduction:- Chern-Simons theories are topological 
quantum field theories that support anyonic particles 
with exotic mutual statistics 1 . In high energy physics 
these theories are encountered in the context of quantum 
anomalies [5J [3] and in the study of gauge theories [I] . 
In the context of condensed matter Chern-Simons the- 
ories emerge as effective theories for the description of 
the fractional quantum Hall liquids 5, 6 J or of the sur- 
face states of three-dimensional topological insulators [7] . 
Double Chern-Simons theories, called BF theories, have 
been recently found application in graphene when it is 
decorated with gauge fields [5]. 

An important high energy physics example that sup- 
ports the Chern-Simons-Maxwell theory is the 2 + 1- 
dimensional massive Thirring model. This model de- 
scribes massive interacting Dirac fermions [9 . It is well 
known that in 1 + 1 dimensions, there exists an exact 
mapping between this massive model and the bosonic 
sine-Gordon model [10] . In the 2 + 1-dimensional case 
the bosonisation gives rise to the Abelian Chern-Simons- 
Maxwell theory in the large fermion mass limit [TT] . 
In this letter we establish a new connection between 
relativistic quantum field theory and condensed matter 
physics. In particular, we derive the 2 + 1-dimensional 
Thirring model from a tight binding model of fermions 
in the following way. It is well known that Dirac 
fermions can faithfully describe the low energy behaviour 
of fermions tunnelling on a honeycomb lattice. An un- 
desirable doubling in these fermionic modes results from 
the lattice nature of the system [12]. Haldane [13] deco- 
rated the honeycomb lattice with next-to-nearest neigh- 
bour tunnelling couplings in such a way that the two 
Dirac modes acquire inequivalent energy gaps. Here 
we employ the adiabatic elimination procedure to freeze 
the dynamics with respect to one of the Dirac modes. 
Subsequently, we introduce interactions between the lat- 
tice fermions and obtain the 2 + 1-dimensional Thirring 



model. While the Haldane model gives rise to the integer 
quantum Hall effect, the interactions introduced in this 
letter are exactly designed to produce fractionalisation 
of charge. Hence, a Chern-Simons theory emerges with 
quasiparticle excitations that are Abelian anyons. This 
theory is accompanied by an additional Maxwell term 
that can be either made negligible or dominant by con- 
trolling the interactions between fermions. It is worth 
noting that, similar to the Haldane model, our model 
breaks time reversal symmetry without a magnetic field. 
Finally, to demonstrate the topological order of the tight 
binding model, we employ the stabilisation of its ground 
state against arbitrary Wilson loop operators. This can 
be shown just by performing local density measurements 
of the lattice fermions. 

From tight-binding to Thirring model:- Let us start 
by describing the fermionic lattice with low energy be- 
haviour given by the 2 + 1-dimcnsional Thirring model. 
Consider the honeycomb lattice, shown in Fig. [IJ with 
fermionic modes placed at each lattice site. Fermions 
tunnel between nearest and next-to-nearest neighbour- 
ing sites. The unit cell of the lattice includes two sites 
that are named b and w. The Hamiltonian of the system 
is given by 

H = -*I>i™j + - E H» 6 i 6 J + e^tflS) 

<ij) «iJ» 

+Uj2f>tbMw i , (1) 
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where t is the nearest neighbour tunnelling coupling and 
if, and t w are the next-to-nearest neighbour tunnelling 
couplings for the b and the w fermions, respectively. Fi- 
nally, U is the interaction coupling that is activated only 
between fermions of the same unit cell. For concreteness, 
we take all couplings, i, %, t w and U to be real and posi- 
tive. A complex phase factor e^ iJ appears explicitly only 
in the next-to-nearest neighbour tunnelling term of the 
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FIG. 1: Left: The honeycomb lattice with fermions tun- 
nelling from one site to the neighbouring (coupling t) as well 
as to the next-to-neighbouring sites (couplings tt and t w ). 
The unit cell is depicted with two sites named b and w. Right: 
The energy dispersion of the t- and it, >u ,-terms of Hamiltonian 
([l]) for momenta that cross trough both Fermi points, P+ and 
P_. The t(,,ii)-term opens up asymmetric energy gaps AE+ 
and AE- to the corresponding Fermi points. We take the 
Fermi energy Ef (depicted with dashed line) between both 
gaps so that only the lower band is completely filled. For 
AE- 3> AB+ we can adiabatically eliminate P_. Here we 
took t=l,t b = 0.02 and t w = 0.1. 



w-particles. Note also the minus phase factor in front of 
the tb couplings. 

For Fermi energies, Ep, close to half filling (see Fig. 
[I]) Hamiltonian ([I]) has the following characteristics. The 
low energy behaviour of the first i-term, H t , is equivalent 
to graphene [14j[T5]. The energy dispersion relation E(p) 
with respect to this term becomes zero for two isolated 
momenta, P± — (0, ±47r/(3v / 3)), called Fermi points. 
Expanding the Hamiltonian around these momenta gives 
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where b(r) and w(r) are the continuous version of the 
fermionic operators and d x , y are partial derivatives in 
the two spatial dimensions. The Hamiltonians Zf t are 
gapless, so they describe massless Dirac fermions. 

The second tf, iU) -term, H tb w , opens an energy gap at 
the Fermi points. We now take the phase acquired 
by w fermions to be = — 27r/3 for the direction 
ni = (3/2, V3/2) (—(j> for the direction — ni) and zero for 
the rest of the directions. Then, close to the two Fermi 
points, i.e. within the low energy approximation, H tb w 
assumes the following forms 

H+ a » -3 J d 2 rt b b{v)^b{v), 

w zj ' d 2 r [t w w{v)^w{v) - i 6 6(r)t&(r)] . (3) 

These Hamiltonians give rise to the energy gaps AE + = 
3tb for P + and AE- = 3(t& + t w ) for P_. Hence, the 
non-zero phase factor <fi allows us to open different gaps 
for the two Fermi points. In particular, we choose 



so the two Fermi points have a large energy difference, 
as shown in Fig. [I] By restricting to low enough energy 
scales, of the order of AE + , the dynamics of P_ will 
be frozen and it can be neglected. To demonstrate this 
consider the ground state, |gs), of the system and two 
excited states, |e+) corresponding to the lowest energy 
excitation at P+ and |e_) corresponding to P_. Next, 
we assign the energy gaps AE + and AE- between each 
of the excited states and the ground state. Assume that 
the system is initially prepared in the ground state |gs). 
Consider a small perturbation in the system that cou- 
ples the ground state to both excited states with equal 
strength of the order of AE + . This perturbation has as 
an effect a negligible population to be transferred to |e_) 
and most of the dynamics to take place only between |gs) 
and |e+). Indeed, by adiabatic elimination we find that 
the maximum population of |e_) state at all times is of 
the order of (fi/ AE-) 2 , which we also verified numeri- 
cally. Hence, we can safely neglect the P- Fermi point 
as long as the perturbations acting on the system satisfy 

n < AE-. 

Finally, the interaction t/-term, Hjj 1 of Hamiltonian 
([T]) is local and acts as a repulsion between the fermions 
in the same unit cell. In the continuous approximation 
it takes the form 
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where for U -C AE- we only consider fermionic modes 
around the P+ Fermi point. Combining all the compo- 
nents together we can write the continuum limit of ([I]), 
up to an overall energy shift, in the following way 
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where ^(r) = (6(r) w(r)) T is the Dirac spinor, are 
the Pauli operators with er = (<x y , a x ), p = (— id x , —id y ), 
jx = Tp^fxip = ip^c y ip, j y = V'TyV' = ip^xip, $ = ip^z 
and jo = ^ip- Moreover, g 2 = ^, c = |t and M = |^|. 
Hence, the nearest neighbour tunnelling coupling corre- 
sponds to the speed of light, the next-to-nearest neigh- 
bour tunnelling coupling gives rise to the mass M of the 
Dirac fermions and the lattice fermion interaction corre- 
sponds directly to the current-current interaction of the 
Thirring model. 

From Thirring model to Chern-Simons theory:- Hamil- 
tonian ^ exactly describes the massive Thirring model 
in 2 + 1 dimensions. We now employ the path integral 
formalism to show the connection of this model to Chern- 
Simons- Maxwell theory [TT]. By applying a Wick rota- 
tion on the temporal coordinate, we can write the cor- 
responding Euclidean partition function of the Thirring 
model as 



tb *C t v 



(4) 



Z Th = / Vil)Vip exp 



^(c^-Mc 2 )f 
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We can introduce a vector field a M through the following 
identity 



cxp (J d?x 9 — 



(8) 



T>a,u exp 



<Tx ( -a M a M + g j% 



so that the exponent of the partition function becomes 
quadratic with respect to the fermionic field. We can 
now integrate out the spinor fields 



r I d' 6 X1p((?+~<J- Mc) 1p 



Dip Dip exp 

exp{ -S eS [a]}, (9) 
and obtain an a^-depcndent effective action given by 



S cS [a] = -clog det (p + - c/L- Mc 



(10) 



Upon applying a Pauli-Villars regularisation [TBI E] to 
the effective action we obtain a parity violating term 



ig 2 Mc 



s ^-^cml dix ^ a ^ + °yMc 



(11) 



which is the Abelian Chern-Simons action up to correc- 
tions of order d/Mc. As we are interested in the be- 
haviour of the ground state of the system, which belongs 
to its low energy sector, the O (d/Mc) and higher order 
terms will have a negligible contribution. Expression ( 11 ) 



comes from one-loop calculations of Feynman diagrams. 
However, the Coleman-Hill theorem [TS] guarantees that 
the Chern-Simons is the dominant term and it receives 
no further concreteness at higher loops. For convenience 
we take Mc/\Mc\ to be positive. 

Next, we introduce an interpolating action Si [a, A] 
given by 



Si[a,A}= (12) 



where A^ is an Abelian gauge field. By integrating the 
partition function of Si [a, A] with respect to A^ or with 
respect to a M it is possible to prove [19] the following 
equivalence between the two different partition functions 



Zt = J Da^DA^e- 8 ^ 



Da^ exp 
VA^exp 



- d 



\8nc 2 / 
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Summing up, through the bosonisation mechanism, we 
have shown that the low energy sector of the 2 + 1- 
dimensional Thirring model is equivalent to a Chern- 
Simons-Maxwell theory. In the standard spacetime with 
Lorentzian signature the topological action appears with 
the coupling 2irc/g 2 . Hence, to enhance the Chern- 
Simons action over the Maxwell term we need to make the 
coupling g 2 small. Alternatively, if we are interested in 
obtaining the electromagnetic action in 2 + 1 dimensions 
then we need to make the coupling g 2 large. Here we are 
interested in the case where the topological action is dom- 
inant. The Chern-Simons term makes the gauge theory 
massive, with a correlation length that decreases propor- 
tionally to g 2 . In particular, the corresponding electric 
and magnetic fields die off exponentially fast away from 
the sources. Nevertheless, the field A v can take non-zero 
values everywhere much like the Aharonov-Bohm effect. 

Measurement of topological order:- Finally, we would 
like to identify the topological order of the tight bind- 
ing model given in ([I]). Initially, we consider the Chern- 
Simons theory. The relevant physical observables should 
be operators that are gauge-invariant as well as metric 
independent. For that we take the Wilson loop opera- 
tors 



W(L) = exp 



A^dx* 



(14) 



where L is an arbitrary link in 2 + 1 spacetime, possibly 
having many strands. It was shown in |20| 12 lj that the 
expectation value of the Wilson loop, (H / (i))csM, can 
be expressed in terms of the linking number known 
also as the Gauss integral, of the link L as (W(L))csm = 
exp ( ± i$ L /(87r)). For a single unknotted loop Lq it 
is = 0, so the expectation value becomes equal to 
(W(L ))csm = 1- H L lies completely on the spatial 
surface of the Chern-Simons theory, i.e. having no time 
component, then this expectation value is evaluated with 
respect to the ground state of the system I^csm) and it 
gives 



(*csm| W(L ) |*csm) = 1. 



(15) 



Hence, the Chern-Simons theory has a ground state that 
is stabilised in terms of Wilson operators of all possi- 
ble loop shapes. For non-trivial ground states or loop 
operators this condition can be satisfied only by states 
that are condensates of all possible loops. Such loop 
condensate states are topologically ordered as they ex- 
hibit non-zero topological entropy [22j ]23j and they have 
non-trivial topological degeneracy when the system is 
wrapped around the torus [25] . These two characteristics 
are the main identification tools of topological order. 
Condition ( 15 ) allows us to determine if the tight- 



(13) 



binding model with a low energy behaviour described by 
the Thirring model is topologically ordered or not. It was 
shown by Fradkin and Schaposnik |11) that the expecta- 
tion value of the Wilson loop can be expressed in terms 
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FIG. 2: A part of the honeycomb lattice with a loop Lq 
depicted that encloses the surface So. The loop encloses a 
number of lattice unit cells, of which the populations b\ 6; 
and tOjWj are measured. When the expectation value of 

exp[i X^igEo (^i^i + w i w i)] tends to a non-zero value for ar- 
bitrarily large areas Eo then the ground state of the system 
is necessarily topologically ordered. 



of fermionic observables of the Thirring model, i.e. 

(W(L))csm 



>Th, 



(16) 



where E is a surface bounded by the loop L. We can 



employ this connection to express condition ( 15 ) of topo- 
logical order in terms of fermionic observables. Consider 
a spatial surface, So > of the Thirring model. The flux of 
the fermionic current through Eo is given in terms of the 
current jq = ip^ip as 



dS^7"V= / dS [&(r) f 6(r) + w(r)t«,(r)] . (17) 



In terms of the tight binding model the flux of the current 
jo through Eo becomes the sum of the fermionic densi- 
ties of both species at the sites enclosed by Eo. Hence, 
the expectation value of the exponential of these popula- 
tions with respect to the ground state of the tight binding 
model, I^tb); is given by 



(^TEslexp 



ies 



I* 



TB, 



= 1. 



(18) 



One can now directly determine if the tight binding 
model is topologically ordered. In Fig. [2] we depict an 
area Eo of the lattice bounded by a loop Lq. Care has 
been taken so that Lq does not cut cells in half as they 
are considered as a single point in space during the con- 



tinuous approximation. Then (18) corresponds to mea- 
suring the populations of b and w fermions, b^b i and 
w i w i , in all sites within the region Eo of the tight binding 
model, constructing their sum and then averaging their 
exponential over different realisations of the lattice ex- 
periment. Note that if the coupling g 2 is large and the 
Maxwell term is dominant over the Chern-Simons action 
then (W(Lo)} w e — K I E °I, w here k is some positive con- 
stant and | Sq I is the area enclosed by the loop Lq [2"5] . 



This quantity decreases exponentially fast as the area of 
the loop is increased. This area law behaviour reveals the 
charge confinement of 2 + 1-dimensional QED !26] and it 
can be directly demonstrated with our scheme. 

Conclusions:- In this letter we have presented a tight 
binding model that gives rise, in the low energy limit, 
to Abclian Chern-Simons theories. We extended a ver- 
sion of Haldane's model, with imbalanced masses of the 
resulting Dirac fermions so that one of them is adia- 
batically eliminated. In this limit interactions between 
fermions exactly reproduce the Thirring model. Upon 
bosonisation the latter model is equivalent to the Abelian 
Chern-Simons theory. A direct method to measure the 
topological order of the system is proposed that requires 
local density measurements of the fermions of the tight 
binding model. These measurements can determine the 
invariance of the ground state under applications of arbi- 
trary Wilson loop operators of the model resulting from 
bosonisation. The generalisation of this procedure to 
non- Abelian Chern-Simons theories [57] is a fundamental 
problem with practical applications to topological quan- 
tum technologies [23] . 

A possible experimental realisation of the tight binding 
model can be given in terms of spin-dependent potentials, 
in the same lines as Refs. [28ll30| . There interspecies 
tunnelling along the honeycomb lattice is activated by 
Raman assisted tunnelling, which can imprint complex 
phase factors as the ones we require here [25]. The in- 
teractions between fermions are restricted only within 
the unit cell, and thus need to be independent of the 
tunnelling couplings. For that one can employ optically 
induced p-wave Feshbach resonance to manipulate the 
collisional couplings U jSU [35] • Finally, the local atom 
density measurements necessary to identify the topolog- 
ical order can be performed with well established tech- 
niques 33J . Relation (18) can then be verified for arbi- 
trary surfaces Eo with geometric characteristics that are 
large compared to the correlation length of the system 

ED- 

Note that an alternative approach to obtain frac- 
tional quantum Hall physics by introducing interactions 
in the Haldane model has been recently presented in [35j , 
though that model is analytically intractable. 
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